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Abstract 

Using a data sample of 57 pb _1 recorded at a centre-of-mass energy of 183 GeV with the Opal 
detector at Lep, 282 W + W~ — ► qqqq and 300 W + W~ — > qqtPi candidate events are used to obtain 
a measurement of the mass of the W boson, J% = 80.39 ± 0.13(stat.) ± 0.05(syst.) GeV, assuming 
the Standard Model relation between M\v and I\y- A second fit provides a direct measure of the 
width of the W boson and gives Tyv = 1.96±0.34(stat.) ±0.20(syst.) GeV. These results are combined 
with previous Opal results to obtain M w = 80.38 ± 0.12(stat.) ± 0.05(syst.) GeV and T w = 1.84 ± 
0.32(stat.) ± 0.20(syst.) GeV. 
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1 Introduction 



Comparison between direct measurements of the mass of the W boson, M\y, and the value determined 
indirectly from precise electroweak results from data taken at y/s ~ Mi and lower energies ]pj provides 
an important test of the Standard Model. In addition, the direct measurement of M\y can be used to 
constrain the mass of the Higgs boson, Me, by comparison with theoretical predictions involving ra- 
diative corrections sensitive to Mh Q. The constraints imposed using Mw, which are presently limited 
by statistical uncertainties, are complementary to those imposed by other electroweak measurements, 
which are largely limited by theoretical uncertainties Q. 

The first direct measurements of the W boson mass were performed at hadron colliders §. During 
1996, the Lep collider at CERN began operating at centre-of-mass energies exceeding the W + W~ 
production threshold (LEP2) thus allowing direct determinations of the mass of the W boson. The 
combination of direct measurements from LEP 2 [5-9] and from hadron colliders presently yields Myv = 



80.41 ± 0.10 GeV [10] . It is expected that LEP2 will ultimately achieve a precision on the W mass of 



approximately 30-40 MeV [jj. 

In 1997, Opal collected 57 pb" 1 of data at a centre-of-mass energy of approximately 183 GeV. 
This paper describes measurements of the W boson mass and width using this data sample. 



2 The OPAL Detector 



The Opal detector includes a 3.7 m diameter tracking volume immersed in a 0.435 T axial magnetic 



field, which yields a transverse^ momentum resolution of & Pxy /p X y ~ y (0.020) 2 + (0.0015 • p xy ) 2 /GeV 2 
and an average angular resolution of about 0.3 mrad in cj) and 1 mrad in 6. The electromagnetic 
calorimeter consists of 11704 lead glass blocks with full azimuthal acceptance in the range |cos#| < 
0.98 and a relative energy resolution of approximately 3% at E ~ 47 GeV, the mean energy of 
electrons from W decays. The magnet return yoke is instrumented with streamer tubes to serve 
as the hadronic calorimeter. Muon chambers surrounding the hadronic calorimeter provide muon 
identification over the range |cos#| < 0.98. Jets are constructed from charged tracks and energy 
deposits in the electromagnetic and hadronic calorimeters using the Durham algorithm [12|. The 
energies of reconstructed jets are calculated using the technique described in [13|. Using e + e~ — > qq 
data taken at = 91 GeV, the jet energy resolution is determined to be approximately &e/E 20% 
with an angular resolution of 20-30 mrad depending on the jet visible energy and polar angle. A more 
detailed description of the Opal detector can be found in 14]. 



3 Data and Monte Carlo Samples 



The integrated luminosity of the data sample, evaluated using small angle Bhabha scattering events 
observed in the silicon tungsten forward calorimeter, is 57.21 ± 0.15(stat.) ± 0.20(syst.) pb _1 . The 
mean centre-of-mass energy, weighted by luminosity, is y/s = 182.68 ± 0.05 GeV | |15| |. 

1 The Opal coordinate system is denned such that the z-axis is parallel to and in the direction of the e~ beam, the 
x-axis lies in the plane of the accelerator pointing towards the centre of the Lep ring, and the y-axis is normal to the 
plane of the accelerator and has its positive direction defined to yield a right-handed coordinate system. The azimuthal 
angle, 0, and the polar angle, 9, are the conventional spherical coordinates. 



3 



3.1 Event Selections 



The event selections are described in detail in | |16| |. The selections are sensitive to the leptonic 
W + W~ — > £ + U£l'~17£i , semi- leptonic W + W~ — > qqtui and hadronic W + W~ — > qqqq final states. 
By construction, the selections are mutually exclusive. The leptonic final state is not used in this 
analysis. 

Semi-leptonic W + W~ — ► qqtvi decays comprise 44% of the total W + W~ cross-section. The 
selection employs three multivariate relative likelihood discriminants, one for each of the W + W~ — ► 
qqeF c , W + W _ — ► qqpV^ and W + W~ — ► qqrV T final states. The W + W~ — > qqez7 e and W + W~ — ► 
qq^v^i channels are characterized by two well-separated hadronic jets, a high- momentum lepton and 
missing momentum due to the prompt neutrino from the leptonic W decay. The W + W~ — > qqtl7 T 
channel is characterized similarly except that the r lepton is identified as an isolated, low-multiplicity 
jet typically consisting of one or three charged tracks. After all cuts, W + W~ — ► qqiPe events are 
selected with an efficiency of 85% and a purity of 90%. The dominant backgrounds are Z°/7 — > qq 
and four-fermion processes such as e + e~ — ► Wel^ c and (Z°/7)*(Z°/7)* — > qq£ + l~. 

Hadronic W + W~ — > qqqq decays comprise 46% of the total W + W~ cross-section and are charac- 
terized by four energetic hadronic jets and little or no missing energy. A loose preselection removes 
approximately 98% of the dominant background process, Z°/7 — ► qq. Following the preselection, a 
multivariate relative likelihood discriminant is employed to select the W + W~ — * qqqq candidates with 
an efficiency of 85% and a purity of 78%. 

After these selections, 361 W + W~ — ► qqiug and 438 W + W~ — ► qqqq candidate events are iden- 



tified, consistent with Standard Model expectations pq| . As discussed in Section |4.1| , additional cuts 
are applied to remove poorly reconstructed events and further reduce backgrounds. 



3.2 Monte Carlo Samples 



A number of Monte Carlo simulation programs are used to provide estimates of efficiencies and purities 
as well as the shapes of the W mass distributions. The majority of the samples are generated at = 
183 GeV assuming Jl% = 80.33 GeV. All Monte Carlo samples include a full simulation of the Opal 
detector [17]. 



The main W + W~ samples are generated using Koralw [18| and include only the CC03 diagrams^. 
The four-fermion backgrounds WeV e and (Z°/7)*(Z°/7)* are simulated using Pythia fll9|| , while grc4f 
[p0[| and Excalibur [21] are used to estimate systematic uncertainties. The background process 
Z°/7 — ► qq is simulated using Pythia, with Herwig [22] used as an alternative to assess possible 
systematic effects. 



4 Measurement of the Mass and Width of the W Boson 



The W boson mass, -Mw, and decay width, Tw, are determined from fits to the reconstructed invariant 
mass spectrum of the selected W + W~ — > qqiVp and W + W _ — > qqqq events. For each selected event, a 
kinematic fit is employed to improve the mass resolution and further reduce background. A reweighting 

2 In this paper, the doubly-resonant W pair production diagrams, i.e. t-channel u c exchange and s-channel Z°/7 
exchange, are referred to as "CC03", following the notation of KM. 
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technique 0] is used to produce Monte Carlo mass spectra corresponding to any given mass and width. 
A binned likelihood fit is used to determine Mw an d T\y by comparing the shape of the reconstructed 
invariant mass distribution from the data to that from reweighted Monte Carlo spectra. 

Two alternative methods are also used to extract Mvv- In the first, an analytic fit to the measured 
mass spectrum uses an unbinned likelihood fit to determine M\y. To describe the signal shape, the fit 
uses a parametrization based on a Breit-Wigner function ||. The second method uses a convolution 
technique similar to that used by the Delphi Collaboration ||. These alternative fits have sensitivities 
similar to that of the reweighting fit and are used as cross-checks. 

4.1 Invariant Mass Reconstruction 

The three methods for extracting Mvv use essentially the same procedures to reconstruct the invariant 
mass of the W candidates. The description provided here applies to the reweighting method. Small 
variations relevant for the alternative analyses are discussed in Section [4.3| . 

For the selected W + W~ — ► qqqq events tracks and clusters are grouped into four jets using the 
Durham algorithm. A kinematic fit is then performed to estimate the reconstructed invariant mass of 
the W candidate. The fit incorporates the constraints of energy and momentum conservation (4C fit) 
yielding two reconstructed masses per event, one for each W boson in the final state. An additional 
constraint can be incorporated by neglecting the finite W width and constraining the masses of the 
two W boson candidates to be equal (5C fit), thus yielding a single reconstructed invariant mass for 
each event. The kinematic fit employs the method of Lagrange multipliers and a ^-minimization 
technique. For extracting il% from the W + W~ — > qqqq candidates, the 5C fit is used to determine 
a reconstructed invariant mass, m rec , its error, a rec , and a \ 2 fit-probability for each event. The 
measured jet momenta with their associated errors and the measured jet masses are used as inputs. 
The use of the measured jet masses, rather than treating the jets as massless, improves the fitted mass 
resolution. Based on Monte Carlo studies, the errors associated with the measured jet momentum are 
parameterized as functions of the visible energy and polar angle of the jet. 

For each W + W~ — > qqqq event three kinematic fits are performed, corresponding to the three pos- 
sible jet pairings. This ambiguity in the choice of the jet pairing leads to a combinatorial background. 
To eliminate poorly reconstructed events and reduce backgrounds, only combinations which yield a 
5C fit with a x 2 fit-probability exceeding 0.01 and m rec > 65 GeV are considered. In addition, combi- 
nations with cj rec < 0.5 GeV are excluded^. A relative likelihood discriminant is employed to choose 
a single combination for each event. The likelihood is constructed for each surviving combination and 
takes as input the following three variables: the difference between the two fitted masses resulting 
from a 4C fit, the sum of the di-jet opening angles and the 5C fit mass. The resulting jet-pairing 
likelihood distribution is shown in Figure |l|. For each event, the combination corresponding to the 
largest jet-pairing likelihood is retained provided it has a likelihood output exceeding 0.18. If the 
combination fails to satisfy this criterion, the event is not used. The cut of 0.18 on the jet-pairing 
likelihood is chosen to optimize the product of efficiency and purity. Approximately 69% of selected 
W + W~ — * qqqq events survive, while 60% of the background is removed. Most of the rejected events 
fail the cut on the \ 2 fit-probability. Monte Carlo studies estimate that in 87% of the surviving 
signal events, the selected combination corresponds to the correct jet pairing and that this fraction 
is independent of J% to within 0.1% over a ±1 GeV range. The number of surviving events in the 
W + W~ — * qqqq channel is given in Table |. 

3 Fits which yield a ICC < 0.5 GeV are excluded because Monte Carlo studies reveal that the reconstructed mass 
resolution of these events is very poor; these events often have a reconstructed invariant mass close to the kinematic 
limit and are assigned an anomalously small fit error. 
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In the selected W + W~ — ► qqel/ c and W + W~ — ► qq^v^ events the non-leptonic part of the event is 
reconstructed as two jets using the Durham algorithm. A kinematic fit is then performed incorporating 
the same five constraints as employed for the W + W _ — * qqqq events. This results in a 2C fit since 
the three-momentum of the neutrino is not known. For the leptons, the inputs to the kinematic fit 
are the lepton energy and direction and their associated errors. The direction is estimated using the 
track associated with the electron or muon candidate. The energy is estimated from the associated 
electromagnetic calorimeter cluster for electrons and from the momentum of the associated track for 
muons. Jets are treated in the manner described above. Only events with a x 2 fit-probability exceeding 
0.001, m rec > 65 GeV and <r rec > 0.5 GeV are retained. These cuts reduce backgrounds by roughly 55% 
and remove poorly reconstructed events. Since the W + W~ — * qqtPe event selections already yield 
low backgrounds, the cut on the x 2 fit-probability is looser than for the W + W _ — > qqqq channel. 
Approximately 90% and 92% of selected W + W~ — ► qqe!7 e and W + W~ — ► qq/^j events, respectively, 
satisfy these additional criteria. Most of the rejected events fail the cut on the x 2 fit-probability. The 
numbers of surviving events are listed in Table p]. 

The Selected W + W~ — > qqrV T events are also reconstructed as two jets using the Durham al- 
gorithm after excluding the tracks and clusters associated with the tau. The invariant mass of the 
jet-jet system, scaled by the ratio of the beam energy to the sum of the jet energies, and its associated 
error are used in determining M\y from this channel. Only events with a reconstructed invariant 
mass greater than 65 GeV and an error on the reconstructed invariant mass greater than 0.5 GeV 
are retained. In addition, to further reduce background, a 1C fit is performed and the resulting x 2 
fit-probability is required to exceed 0.001. The fit incorporates the same five constraints as employed 
for the W + W~ — > qqqq channel, assumes that the r lepton direction is given by the direction of the 
visible decay products associated with the tau and estimates the total energy of the tau using energy 
and momentum constraints. Approximately 78% of selected W + W _ — > qqr!7 T events satisfy these 
additional criteria while 65% of the background is removed. Most of the rejected events fail the cut 
on the x 2 fit-probability. The number of surviving events is listed in Table [|. 

The full width at half maximum of the residual of the reconstructed invariant mass per event is 
used as an estimate of the average m Tec resolution and is calculated using Monte Carlo events in which 
less than 100 MeV of energy is radiated into initial state photons. For W + W~ — ► qqqq events this 
resolution is 1.7 GeV for fits corresponding to the correct jet pairing. For W + W~ — ► qqtPe events the 
average m rec resolution per event is 2.4, 2.8 and 3.4 GeV in the W + W~ — > qqeUc, W + W~ — > qq/^V 
and W + W~ — * qqru T channels, respectively. 

4.2 Extraction of the W Mass and Width 

The W boson mass and width are extracted by directly comparing the reconstructed mass distribution 
of the data to mass spectra obtained from fully simulated Monte Carlo events corresponding to various 
values of JV% and I\v. A likelihood fit is used to extract Myv and Tyv by determining which Monte 
Carlo spectrum best describes the data. The Monte Carlo spectra for arbitrary values of Mw an d T\y 
are obtained using the Monte Carlo reweighting technique described in fj. 

The mass spectra for background events are taken from Monte Carlo and are assumed to be 
independent of Jl% and Tw- The background reconstructed mass distributions are normalised to the 
expected number of background events. The reweighted signal spectra are then normalised such that 
the total number of signal plus background events corresponds to the observed number of events. This 
is done separately for the W + W~ — ► qqqq, W + W~ — > qqeI7 e , W + W~ — ► qq/tt 7 ^ and W + W~ — > qqrz7 T 
channels. In addition, the W + W~ — * qqive channels are divided into four subsamples according to the 
error on the reconstructed invariant mass. These subsamples are treated independently within each 
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channel. This division gives a larger weight to events with reconstructed masses which are known 
with better precision (i.e. small cr rec ) relative to events with poorly determined m rec and reduces 
the expected statistical uncertainty on the fitted W mass by approximately 7% in the W + W~ — * 
qqtug channels. In the W + W~ — * qqqq channel, the width of the reconstructed mass distribution is 
dominated by the intrinsic width of the W so that a similar subdivision does not improve the Mw 
sensitivity in this channel and therefore is not implemented. The <7 rec distribution is shown in Figure |2| 
for the W + W~ — > qqqq and W + W~ — > qqiUg channels separately. 

A binned log-likelihood fit to the m rec distributions of the data is performed in the range m rec > 
65 GeV. The log-likelihood function is defined as 



where n« is the number of observed events in the ith bin, fj, is the expected background fraction 
of the sample using the normalization procedure described above, ^(Mwdw) is the probability of 
observing a signal event in the ith bin assuming a W boson mass and width of Myy and T\y and V\ 
is the analogous probability for the background, which is assumed to be independent of the W mass 
and width. Both and V\ are estimated using the relevant Monte Carlo spectrum. The 

log-likelihood curve is determined separately for each channel. For the W + W~ — > qqi^i channels, the 
results are obtained by adding the log-likelihood curves separately determined from each subsample. 

Two types of fit are performed. In the one-parameter fit, Tw is constrained to its Standard Model 
relation to the W mass H and only Jl% is determined. The results of this fit for each channel are 
given in Table || and displayed in Figure ||. The combined result is discussed in Section ^. In the two- 
parameter fit, both Mw and Tw are determined simultaneously. The likelihood contours for this fit 
are displayed in Figure || including statistical errors only. The systematic uncertainties are discussed 



One advantage of the reweighting method is that the fitted parameters should be unbiased since 
any offsets introduced in the analysis are implicitly accounted for in the Monte Carlo spectra used in 
the reweighting procedure. This is verified using several Monte Carlo samples generated at various M\y 
and In addition, tests using a large ensemble of Monte Carlo subsamples, each corresponding to 
57 pb _1 and including background contributions, are used to verify for each channel separately and for 
all channels combined, that the measured fit errors accurately reflect the RMS spread of the residual 
distribution for both the J% and T\v fits. 

4.3 Alternative Fit Methods 
4.3.1 Breit-Wigner Fit 

The Breit-Wigner method is analogous to that described in ||. It employs an unbinned maximum- 
likelihood fit to the reconstructed mass spectrum using an analytic Breit-Wigner function to describe 
the signal. Due to initial-state radiation, the reconstructed mass spectrum is significantly asymmetric 
for data taken at y/s ~ 183 GeV. It is found that a relativistic Breit-Wigner function, with different 
widths above and below the peak, gives a satisfactory description of the m rec lineshape. The function 
is given by 




i=l 



in Section [|. 



S (m re c) = -4 





'rec 



m l) 2 + m rec T - 



/+ 
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where r+(— ) is the width assumed for all m TCC above (below) the peak, tuq. The widths are fixed to 
values determined from fits to W + W~ signal Monte Carlo samples and are found to be independent 
of M\v over the range relevant for this analysis. The shapes of the background distributions and the 
background fraction are also determined from Monte Carlo. The background fraction is held constant 
in the fit. The fit is performed over the range 70 < m Tec < 88 GeV. 

For the W + W~ — ► qqqq events, the likelihood method for choosing which jet pairing to use is 
found to distort the m Tec distribution so that it is inadequately described by a Breit-Wigner function. 
Therefore, the following procedure is employed ||. The reconstructed invariant mass of the combi- 
nation with the largest fit probability, Pi, is used if -Pi > 0.01. The reconstructed invariant mass of 
the combination with the second largest fit probability, P2, is also used if it satisfies P% > 0.01 and 
P2 > Pl/3. The selected reconstructed masses enter the same distribution with unit weight. Monte 
Carlo studies estimate that the correct combination is among those chosen in 90% of the surviving 
W + W~ — > qqqq events. 

In contrast to the procedure employed for the reweighting method, the W + W~ — > qqtPg events are 
not divided into subsamples according to a rec because the subsamples exhibit a distorted reconstructed 
mass distribution which is poorly described by a Breit-Wigner function. 

The fitted mass, mo, must be corrected for offsets not accounted for in the fit, e.g. from initial- 
state radiation and event selection. A correction is determined using fully simulated Monte Carlo 
samples generated at various known M\y and yfs with the expected background contributions included 
and is found to depend linearly on both Myy and yfs. The results from the W + W~ — ► qqqq and 
W + W~ — ► qqiVi channels, after correction, are given in Table ||. 

4.3.2 Convolution Fit 

The convolution method is similar to that employed by the Delphi Collaboration [||. The method 
attempts to exploit all available information by constructing a likelihood curve for each selected event. 
The likelihood is calculated using the functional expression 

where p s is the probability of the candidate event being a real signal event and V s is the probability 
density function for the signal, 

V S (M W , m rec ) = B W(M W , m) tg> R(m, m rec ), 

where R(m,m rec ) is the resolution function estimated from the Monte Carlo and BW(il%,m) is a 
relativistic Breit-Wigner function with Tw fixed to its Standard Model relation to Myv- The exact 
expression for the likelihood is channel dependent. For example, in the W + W~ — > qqqq channel 
the added complication of combinatorial background requires a sum over the three jet pairings. The 
log-likelihood curves from each selected event are summed to yield a single curve from which a fitted 
mass is determined. This fitted mass is corrected for background effects and for offsets not accounted 
for in the fit in the manner described for the Breit-Wigner fit. The results from the W + W~ — > qqqq 
and W + W~ — > qqtPi channels, after all corrections, are given in Table |. 

5 Systematic Uncertainties 

The systematic uncertainties are estimated as described below and summarised in Table H 
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Beam Energy: 



The average Lep beam energy is known to a precision of ±25 MeV [15|. The effect of this 
uncertainty on the measured M\y is determined from fits to large Monte Carlo samples for 
which the analysis is repeated assuming = ± 50 MeV. The observed shifts in the fitted 
M\v are used to estimate the associated systematic uncertainty. The spread in Lep beam energy 



of 152 ± 8 MeV [15| has a negligible effect on both the mass and width determinations. 



Initial State Radiation: 

The systematic error associated with uncertainties in the modelling of initial state radiation is 
estimated by comparing fully simulated Koralw Monte Carlo W + W~ events generated using 
a leading logarithm 0(a) treatment of initial state radiation to the standard W + W~ sample 
which includes a next-to-leading-log 0(a 2 ) treatment. No significant difference is observed and 
the statistical uncertainty of the comparison is taken as the associated systematic error. 

Hadronization: 

The scale of hadronization effects is studied by comparing the fit results of a single W + W~ 
sample generated once using Pythia and again using Herwig as the hadronization model. Both 
samples contain the same W + W~ final states and differ only in their hadronization modelling. No 
significant differences are observed and the statistical uncertainty of the comparison is taken as 
the associated systematic error. As a cross-check, W + W _ samples are generated with variations 
of the Jetset fragmentation parameters a q , b, Aqcd and Qq, of one standard deviation about 



their tuned values |23|]. The fit results are compared and yield no statistically significant effects. 



Four-fermion Effects: 

The Monte Carlo samples used to estimate background contributions in the reweighting pro- 
cedure do not include a complete set of four-fermion diagrams and neglect interference effects 
between W + W~ diagrams and other four-fermion processes. In order to test the sensitivity of 
the results to these effects, the fit results of a sample generated including the full set of interfering 
four-fermion diagrams are compared to those of a sample restricted to the CC03 W + W~ dia- 
grams alone. The comparison is performed using both the grc4f and the Excalibur generators. 
In neither case is a significant difference observed. The larger of the two statistical uncertainties 
is assigned as the associated systematic error. 

Detector Effects: 

The effects of detector mis-calibrations and deficiencies in the Monte Carlo simulation of the data 
are investigated by varying the jet and lepton energy scales and the errors input to the kinematic 
fit over reasonable ranges. The ranges are dependent on polar angle and are determined from 
detailed comparisons between 1997 data and Monte Carlo using approximately 2.1 pb" 1 of data 
collected at Mz and e + e~ — > e + e~ events recorded at y/s ~ 183 GeV. Of particular 

importance for the W + W _ — * qqft^ channel are the lepton energy scales, which are varied by 
±0.5%. For the W + W~ — > qqqq channels the most important variations are for the errors 
associated with the jet angles, which are varied by ±7% for errors in cos 9 and ±3% for errors 
in <f). The jet energy scale is varied by ±1-2% over most of |cos0|. For each variation a large 
Monte Carlo sample is refitted and the resulting shifts in the fitted Mw are added in quadrature 
to yield an estimate of the associated systematic error. 

Fit Procedure: 

The reweighting procedure accounts for the fact that the Monte Carlo samples are generated at 
a centre-of-mass energy different from that of the data ||. Monte Carlo samples generated at 
yfs = 182 and 184 GeV are used to test this procedure and assign a systematic error. 

The three methods used to measure M\\r each utilise the data and Monte Carlo in different 
ways. A comparison of the fitted results is used to test for residual biases in the reweighting 
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fit. The three methods are compared using an ensemble of Monte Carlo subsamples, each 
corresponding to 57 pb _1 and including background. For each subsample the difference in the 
fitted Jl% determined using the reweighting method and that determined using each of the 
other methods is calculated for the W + W~ — > qqqq and W + W~ — ► qqfi^ channels separately 
and for the combined sample. The distributions of these differences are approximately Gaussian 
and have means consistent with zero and RMS values of approximately 90 MeV, 95 MeV, and 
75 MeV (100 MeV, 75 MeV, and 60 MeV) when comparing the reweighting fits to the Breit- 
Wigner (convolution) fits in the W + W~ — > qqqq, W + W~ — > qqiug and the combined samples, 
respectively. Since these alternative analyses yield results consistent with those obtained using 
the default reweighting analysis, no additional systematic is assigned based on these comparisons. 

Background Treatment: 

Uncertainties associated with both the normalization and shape of the background distributions 
are investigated. The background normalization is varied by one standard deviation of its asso- 
ciated uncertainties as evaluated in and the data are refitted. As an estimate of the errors 
associated with the uncertainties in the shape of the background distributions, a variety of sub- 
stitutions are made. For all channels, Herwig replaced Pythia for the hadronization model 
for the Z°/7 — ► qq background. In addition, for the W + W~ — > qqqq channel, data taken at 
y/s Mz, scaled by (183 GeV/Mz), are also substituted for the Z°/7 — > qq background. For 
each substitution the data are refitted. The quadrature sum of the shift in fitted mass observed 
when changing the normalization and the largest of the observed shifts in the fitted mass from 
the various substitutions in each channel is assigned as a systematic error. Using Pythia sam- 
ples of e + e~ — > WeU e events generated using various M\y, it is verified that the uncertainty on 
M\v has a negligible effect on the background distributions. 

Monte Carlo Statistics: 

The finite statistics of the Monte Carlo samples used in the reweighting procedure contribute a 
systematic uncertainty of ±15 MeV to the W mass determined separately in the W + W~ — * qqqq 
and W + W~ — ► qqiPe channels, and ±40 MeV to the W width determined from the combined 
sample. 

Colour-Reconnection Effects and Bose-Einstein Correlations: 



As discussed in [11] and [|2J] and references therein, a significant bias to the apparent W mass 
measured in the W + W~ — ► qqqq channel could arise from the effects of colour-reconnection 
and /or Bose-Einstein correlations between the decay products of the W + and W - . These effects 
are investigated separately. 

Using currently available W + W~ data, it is not possible to discern whether or not Bose-Einstein 
correlations are present between hadrons originating from different W decays [25-26]. To in- 
vestigate possible systematic biases a Monte Carlo sample is generated including Bose-Einstein 



correlations using the Pythia Monte Carlo generator implemented as described in [27]. The 
fit result in the W + W~ — > qqqq channel from this sample is compared to a fit from a PYTHIA 
sample which excludes Bose-Einstein correlations. No significant bias is observed for the fitted 
parameters and the statistical precision of the comparison is taken as the associated systematic 
error of ±32 MeV in the fitted mass determined from the W + W~ — * qqqq channel. The uncer- 
tainty on the fitted width in W + W~ — ► qqqq channel propagates to an uncertainty of ±55 MeV 
on the fitted width determined from the combined sample. 

To investigate the systematic biases originating from colour-reconnection effects, several models 
are studied using the Pythia and Ariadne p3| Monte Carlo generators. As discussed in the 



accompanying paper [29], comparing various event shapes offers a means of discriminating be- 
tween the models and, when comparing with data, a means of testing each model independently 
of the M\v measurement. Based on these studies, and comparing fitted masses between Monte 
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Carlo samples including and excluding colour-reconnection effects, a systematic uncertainty of 
±49 MeV is assigned to the mass determined from the W + W~ — ► qqqq channel |29|. The uncer- 
tainty on the fitted width in W + W~ — > qqqq channel propagates to an uncertainty of ±109 MeV 
on the fitted width determined from the combined sample. 



The contributions from each of the above sources are added in quadrature to yield the total systematic 
uncertainty. For the alternative analyses, the systematics are estimated similarly and yield comparable 
results. 



6 Results 



For the reweighting method described in Section 4.2, the results of a simultaneous fit to Mw and T\y 



from the combined W + W — > qqqq and W + W — ► qqiui event samples are 

M w = 80.39 ± 0.13 ± 0.06 GeV, 
T w = 1.96 ±0.34 ±0.20 GeV, 

where the uncertainties are statistical and systematic, respectively. The correlation coefficient between 
Mw and Tw is 0.13. For this fit, the central values are determined by adding the log-likelihood curves 
from the W + W~ — > qqqq and W + W~ — > qqi^i channels. The systematic uncertainties are also 
estimated by summing the log-likelihood curves from each channel. 

A one-parameter fit for the mass is performed by constraining the width using the Standard 
Model relation to give M w = 80.53 ± 0.23(stat) ± 0.09(syst) GeV in the W+W" -> qqqq channel, 
and M w = 80.33 ± 0.17(stat) ± 0.06(syst) GeV in the W+W" -> qq£v e channel. The combined 
result is determined taking into account the correlated systematics between the W + W~ — > qqqq and 
W + W~ — > qqiPi channels and gives 

M w = 80.39 ±0.13 ±0.05 GeV, 

where the uncertainties are statistical and systematic, respectively. For the combination, the W + W~ — > 
qqqq channel carries a weight of 0.34. The combined W + W~ — > qqqq and W + W~ — ► qqd>t re- 
sults from the alternative analyses, after all corrections, are for the Breit-Wigner fit, Mw = 80.37 ± 
0.15(stat) ± 0.05(syst) GeV, and for the convolution fit, M w = 80.30 ± O.H(stat) ± 0.06(syst) GeV. 
As discussed in Section [|, these results are statistically consistent with those obtained using the 
reweighting fit. 

The difference between the fitted Mw in the W + W~ — * qqqq and W + W~ — * qqjtvi channels is 
AM W = (Mffi® - M^ Wl ) = 0.20 ± 0.28 ± 0.07 GeV, where the uncertainties are statistical and 
systematic (excluding contributions from colour-reconnection/Bose-Einstein effects), respectively. A 
significant non-zero value for AMw could indicate that colour-reconnection/Bose-Einstein effects are 
biasing the Mw determined from W + W~ — > qqqq events. 



6.1 Combination with Previous Data 



The measurements of Mw from direct reconstruction at y/s ~ 172 GeV ^ and \/s ~ 183 GeV 
are combined with the Mw measurement from the W + W~ production cross-section at threshold, 
\fs ~ 161 GeV [||. The combination is made assuming that the mass measurements from direct 
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reconstruction and from the threshold cross-section are uncorrelated, apart from the uncertainty as- 
sociated with the Lep beam energy, which is taken to be fully correlated. The direct reconstruction 
measurements are combined accounting for correlated systematics. The combined result is 

M w = 80.38 ± 0.12 ± 0.04 ± 0.02 ± 0.02 GeV, 

where the uncertainties are statistical, systematic, colour-reconnection/Bose-Einstein and beam en- 
ergy, respectively. 

The measurements of T\y from direct reconstruction at y/s sa 172 GeV || and yfs as 183 GeV are 
combined taking into account the correlated systematics to obtain 

T w = 1-84 ± 0.32 ± 0.15 ± 0.12 ± 0.01 GeV, 

where the uncertainties are statistical, systematic, colour-reconnection/Bose-Einstein and beam en- 
ergy, respectively. Monte Carlo studies reveal that the measured statistical error is correlated with the 
measured width. To avoid biasing the combination, the separate width measurements are weighted 
using the expected statistical error, determined from an ensemble of many Monte Carlo experiments 
assuming a W width of 2.093 GeV. 

The measurements of AMw from direct reconstruction at y/s ~ 172 GeV and y/s ~ 183 GeV are 
combined taking into account the correlated systematics to obtain 

AM W = 0.08 ± 0.26 ± 0.07 GeV, 

where the uncertainties are statistical and systematic, respectively. Note that no systematic error is 
included for uncertainties in the modelling of colour-reconnection/Bose-Einstein effects and that the 
uncertainty in the Lep beam energy does not contribute a systematic error to this quantity. 



7 Summary 

Using the 57 pb _1 of data recorded by the Opal detector at a mean centre-of-mass energy of approx- 
imately 183 GeV, a total of 582 W + W~ — > qqqq and W + W~ — > qqiUi candidate events are used in 
a fit constraining Tw to its Standard Model relation with Jl% to obtain a direct measurement of the 
W boson mass, Mw = 80.39 ± 0.13(stat.) ± 0.05(syst.) GeV, while a second fit is used to directly 
determine the width of the W boson, T w = 1-96 ± 0.34(stat.) ± 0.20(syst.) GeV. 

The results described in this paper are combined with the previous Opal results from data recorded 
at y/s ~ 172 GeV and y/s ~ 161 GeV. From this combined data sample the W boson mass is 
determined to be 

M w = 80.38 ± 0.12 ± 0.05 GeV. 

The result for the W boson width is combined with the previous Opal result from data recorded at 
y/s Pa 172 GeV to obtain 

T w = 1-84 ± 0.32 ± 0.20 GeV. 
The uncertainties are statistical and systematic, respectively. 
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Channel 


Observed 


Expected 


Purity 


W+W" -» 


qqqq 


282 


278.1 


86% 


w+w-^ 


qqeI7 e 


119 


113.0 


97% 


w+w-^ 


qq/^V 


107 


114.4 


99% 


w+w-^ 


qqrV T 


74 


82.6 


95% 


Combined 




582 


588.0 


92% 



Table 1: Numbers of events used in the W mass and width determination for each channel and all 



channels combined. Only events surviving the cuts described in Section 4.1 are included. The number 
of expected events and corresponding purities are estimated assuming the world average A% and have 
relative uncertainties of approximately 3%, dominated by the uncertainty in the CC03 production 
cross-section. 



Channel 


Measured M\y Expected error 


W + W _ — >■ qqeI7 e 
W + W~ — > qqfJiPfj, 
W + W~ — > qqrv T 


80.23 ±0.24 0.31 
80.42 ±0.27 0.30 
80.40 ± 0.40 0.45 


W+W~ qqtVz 
W + W~ — > qqqq 


80.33 ±0.17 0.19 
80.53 ±0.23 0.20 



Table 2: Results using the reweighting method for the fit constraining Tyj to its Standard model 
relation with M\y from 57 pb _1 of data taken at yfs « 183 GeV for each of the channels separately 
and for the combined W + W~ — > qqi^e channel. The expected errors are estimated using an ensem- 
ble of Monte Carlo subsamples, each corresponding to 57 pb _1 of data and including background 
contributions. The errors are statistical only. (All quantities are in GeV.) 



Channel 


Breit-Wigner fit 
Measured Mw 


Convolution fit 
Measured M\y 


W+W" -» qqU7 e 
W + W" — ► qqqq 


80.27 ±0.19 
80.52 ±0.24 


80.24 ±0.19 
80.38 ±0.21 



Table 3: Fit results using the alternative analyses and 57 pb _1 of data taken at y/s ~ 183 GeV for 
the W + W~ — > qqiUi and W + W _ — > qqqq channels separately. The expected statistical errors are 
very similar to those of the reweighting method given in Table |2[ The errors are statistical only. (All 
quantities are in GeV.) 
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Systematic errors 
(MeV) 


qqqq 


M w 


comb. 


r w 


AM W 


Beam Energy 




22 


22 


22 


5 





Initial State Radiation 




10 


10 


10 


15 


10 


Hadronization 




21 


21 


16 


52 


15 


Four-fermion 




30 


28 


21 


52 


40 


Detector Effects 




38 


31 


26 


90 


46 


Fit Procedure 




15 


15 


15 


45 


20 


Background 




25 


10 


10 


79 


25 


MC statistics 




15 


15 


11 


40 


21 


Sub-total 




67 


58 


49 


154 


74 


Bose-Einstein Correlations 




32 





11 


55 




Colour Reconnection 




49 





16 


109 




Total systematic error 


89 


58 


53 


196 


74 



Table 4: Summary of the systematic uncertainties for the fit results. For the fits to determine Mw, 
T-yv is constrained to its Standard Model relation. The uncertainties are given separately for fits to 
the W + W~ — > qqqq, W + W~ — > qqii^e and the combined samples. For Tw the uncertainties are given 
only for the fit to the combined sample. The quantity AMw = ^M^™ — M^ 57 ^ uses the Standard 

Model constrained fits to the W + W~ — > qqqq and W + W~ — > qq£ve channels separately. Since the 
primary interest in the quantity AMw is as a test of possible colour-reconnection/Bose-Einstein effects 
in the W + W~ — > qqqq channel, no systematic error is assigned for uncertainties associated with the 
modelling of these effects. The inter-channel correlations are taken into account for the combined Mw 
and Tw fits and for the AMw determination. 
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Figure 1: The jet-pairing likelihood distribution for selected W + W~ — > qqqq events. For each event, 
only the likelihood output for the combination yielding the maximum likelihood is plotted. Events to 
the right of the arrow are retained for the mass analysis. The points correspond to the Opal data and 
the open histogram to the Monte Carlo prediction. The contribution from the non-WW background is 
shown as the cross-hatched histogram and the addition of the combinatorial background is indicated 
by the singly-hatched histogram. 
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Figure 2: The distribution of the error on the reconstructed invariant mass, o" rcc , for W + W~ — > qqqq 
events surviving the jet-pairing likelihood cut and used in the mass analysis (top), and for all selected 
W + W _ — ► qqfi^ events (bottom). The W + W~ — > qq^£ events to the right of the arrow are retained 
for the mass analysis. This error is used to divide the W + W~ — ► qqiUe events into four subsamples: 
0.5 < cr rcc < 1.5 GeV, 1.5 < cr rcc < 2.5 GeV, 2.5 < cr rcc < 3.5 GeV, a TCC > 3.5 GeV, which are 
chosen to minimize the expected statistical error on Jl%. The points correspond to the Opal data 
and the open histogram to the Monte Carlo prediction. The background contribution is indicated 
by the cross-hatched histogram. In the W + W~ — > qqqq channel, the addition of the combinatorial 
background is indicated by the singly-hatched histogram. 
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Figure 3: The reconstructed invariant mass distribution for the W + W~ —> qqqq, W + W~ — »■ qqeI7 e , 
W + W~ — > qq/il 7 ^ and W + W~ — > qqr!7 r samples. The points correspond to the Opal data and 
the open histogram to the reweighted Monte Carlo spectrum corresponding to the fitted mass. The 
background contribution is indicated by the cross-hatched histogram. 
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Figure 4: The 39% and 86% contour levels of the two-parameter fit using the reweighting method. 
The projections of these contours onto the axes give the one and two standard deviation statistical 
uncertainties. The one standard deviation region, including only the statistical error, of the fit to M\y 
only is given by the shaded band. This fit is constrained to the solid line, which gives the dependence 
of the width on the mass according to the Standard Model. The Standard Model prediction for T\y, 
assuming the world average M\y , is shown as an "X" . 
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